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CN ' Abstract 
> 

lO . The search for D^D^ mixing may carry a large discovery potential for new 

physics since the D^D'^ mixing rate is expected to be small in the Standard 
Model. The past decade has seen significant experimental progress in sensitivity. 
I This paper discusses the techniques, current experimental status, and future 

l/^ ■ prospects for the mixing search. Some new ideas, applicable to future mixing 

0^ I searches, are introduced. In this paper, the importance of separately measuring 

i the decay rate difference and the mass difference of the two CP eigenstates 

i' (in order to observe New Physics) has been emphasized, since the theoretical 

calculations for long distance effects are still plagued by large uncertainties. 
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1 Introduction 



Particle-antiparticle mixing has always been of fundamental importance in testing the 
■ Standard Model and constraining new physics. This is because mixing is responsible for the 

small mass differences between the mass eigenstates of neutral mesons. Being a flavor chang- 
ing neutral current (FCNC) process, it often involves heavy quarks in loops. Such higher 
order processes are of great interest since the amplitudes are sensitive to any weakly-coupling 
quark flavor running around the loop. Historically, K^K^ mixing is the rare (FCNC) pro- 
cess that has been experimentally examined in the greatest detail. It has been amply 
demonstrated that in spite of many inherent uncertainties of strong interaction physics, the 
Standard Model predicts the correct phenomenology of the K^K^ mixing. In fact, based on 
the calculation of the Kl - Ks mass difference, Gaillard and Lee were able to estimate 
the value of the charm quark mass before the discovery of charm. Moreover, B^B^ mixing 
gave the first indication of a large top quark mass. 

Although D^D'^ mixing is very similar to K^K^ and B^B'^ mixing, as all are FCNC 
processes, there are significant differences which make D'^D'^ mixing a possible unique place 

Presented at the r-charm Factory Workshop, Argonne National Laboratory, June 20-23, 
1995. 
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to explore new physics. Roughly speaking, in the case of K and B FCNC processes, the 
appearance of the top quark in the internal loop with rat > Mw >> rricrnu removes the 
GIM § suppression, making K and B decays a nice place to test FCNC transitions and 
to study the physics of the top. In the case of D FCNC processes, the FCNC are much 
stronger suppressed because the down- type quarks {d, s and b) with m^, mg, mi, << Mw 
enter the internal loops and the GIM mechanism is much more effective [Q. Therefore 

the D"D mixing rate is expected to be small in the Standard Model, which means the 
mixing search may carry a large potential for discovery of new physics. There are many 
extensions of the Standard Model which allow D'^D mixing (the mass difference between 
the two CP eigenstates) to be significantly larger than the Standard Model prediction (for 
example, see p to pq|). Recent reviews on FCNC processes in D decays can be found 



elsewhere ||16|, |17| |18|, 19 . In general, there could be a large enhancement of the one-loop 
induced FCNC processes in D decays with no constraint from limits on FCNC processes 
in the K and B systems. Roughly speaking, this is because the couplings of FCNC to 
up-type quarks (u,c,i) could be completely different from those to down-type quarks {d,s,b). 
Thus one gains independent pieces of information when searching for FCNC in D decays, 
compared to what is learned searching for FCNC in K and B decays. 

One can characterize D^D^ mixing in terms of two dimensionless variables: x = Sm/j^ 
and y ~ 7-/7+, where the quantities 7± and Sm are defined by j± = (71 ±72)72 and 
Sm = 7712 — mi with mi,"fi (i = 1,2) being the masses and decay rates of the two CP 
(even and odd) eigenstates. Assuming a small mixing, namely, (5to,7_ <C 7+ or x,^/ ^ 1, 
we have Rmixing = {x'^ + 2/^)/2- Mixing can be caused either by x 7^ (meaning that 
mixing is genuinely caused by the — transition) or by ?/ 7^ (meaning mixing is 
caused by the fact that the fast decaying component quickly disappears, leaving the slow 
decaying component which is a mixture of D'^ and i?"). Theoretical calculations of D'^D^ 
mixing in the Standard Model are plagued by large uncertainties. While short distance 
effects from box diagrams are known |^ to give a negligible contribution (^ 10^^°), the 
long distance effects from second-order weak interactions with mesonic intermediate states 
may give a much larger contribution. Estimates of Rmixing from long distance effects range 
from 10~^ to 10"'^ It has recently been argued by Georgi and others that the long 

distance contributions are smaller than previously estimated, implying that cancellations 
occur between contributions from different classes of intermediate mesonic states 1 22| . While 
many people now believe that within the Standard Model Rmixing < 10^^ |l^, Tsj, others 
think Rmixing could be much larger say lO"** |Q (meaning both x and y are above 

10^''). For example, Bigi [ p3[ pointed out that observing a non- vanishing value for Rmixing 
between 10""* and 10^'^ would at present not constitute irrefutable evidence for New Physics, 
considering the large uncertainties in the long distance calculations. While there is some 
hope that the uncertainties can be reduced in the future, as pointed out by Bigi |23 , partly 
through theoretical efforts and partly through more precise and comprehensive data (since 
a more reliable estimate can be obtained from a dispersion relation involving the measured 
branching ratios for the channels common to and D'^ decays), one recent paper claims 



that the hope is rather remote |24 . Speculations abound, but (fortunately) physics is an 
experimental science, and only with solid experimental evidence will we be able to properly 
address these problems. As experimentalist, I think the best way is to measure x and y 



separately, as suggested in |26[ g7|. As will be discussed, this is experimentally possible. If 
we can measure Rmixing as well as y, then we can in effect measure x. Within the Standard 
Model, X and y are expected to be at the same level, although we do not know exactly at 
what level as theoretical calculations for long distance effects (which contribute to both x 
and y) are still plagued by large uncertainties. We expect New Physics does not affect the 
decays in a significant way thus does not contribute to y, but only to x. The point I am 
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trying to make here is that the long distance contribution can be measured, even if it cannot 
be calculated in a reliable way; that is, by measuring y directly. If we can experimentally 
confirm that indeed x » y, then we can claim New Physics, regardless of what theoretical 
calculations for long distance effects are. Otherwise, if it turns out x ^ y, then mostly likely 
we are seeing the Standard Model Physics. Therefore, it is crucial to measure y in order to 
understand the size of x within the Standard Model. This is one of the major points I have 
been trying to make in the past ^ and in this paper. 

Motivated by the experience with K'^K^ system, experimenters have been searching 
for D^D*^ mixing since shortly after the discovery of D° meson at SPEAR in 1976, in 
either hadronic decays D° K~^'k~{X) j28j, or semileptonic decays D'^ 

X^l^v. The past decade has seen significant experimental progress in sensitivity (from 
20% to 0.37% H to @), as can be seen in Figure 1. The search for D°Z)0 mixing has 
a long and interesting history (see Figure [^). In the first few years, people searched for 
— > X+tt" assuming that it would be due to mixing only. Normally, Z?" decays by 
Cabibbo favored decay — > K~tt~^ and 5" — > _ftr+7r~. A signal for if+Tr" could 

indicate mixing of Z)°. But it could also indicate a different decay channel, namely, 

Doubly Cabibbo Suppressed Decay(DCSD) _D" K^n^ , which is suppressed with respect 
to the Cabibbo favored decay by a factor of tan^Oc ^ 0.3% where 9c is the Cabibbo angle. 
As will be discussed, around 1985 there were hints of — > K^-k~-k^ observation, which 
could be due to DCSD or mixing. The popular interpretation neglected the possible DCSD 
contribution, giving the impression that D^D'^ mixing rate Rmixing could be of order 0(1%). 
This engendered much theoretical work to accommodate the possibly large mixing rate. At 
that time, the "theoretical prejudice" was that long-distance contributions dominated and 
would give a large mixing rate on the order of 1% level. Later on, fixed target experiments 
published limits which were not much larger than the nai've quark model DCSD rate. In 
light of these results, the commonly held impression was then that DCSD was much larger 
than mixing, and that exploring mixing by means of hadronic decays had been almost 
exhausted as a technique since the "annoying DCSD background" would inherently limit 
ones ability to observe the interesting physics - D^D^ mixing. It was believed by many that 
the signature of mixing appears only at longer decay times; therefore, it will suffer from 
DCSD fluctuation, and destructive interference could wipe out the signature of mixing. Since 
semileptonic decays are not subject to this "annoying background" , the general consensus 
was that semileptonic decays were a better avenue to explore D^D*^ mixing. 

However, as will be discussed in more detail later, the commonly believed "annoying 
DCSD background" does not necessary inherently limit the hadronic method as the poten- 
tially small mixing signature could show up in the interference term |26j. Moreover, the 
possible differences between the resonant substructure in many DCSD and mixing decay 
modes could, in principle, be used to distinguish between DCSD and mixing candidates 
experimentally (the importance of the mixing-DCSD interference effect will be more 
clear here). Our ability to observe the signature of a potentially small mixing signal de- 
pends on the number of D*^ K^Tr~ (X) events we will have. This means observing 
— > K^Tr~ (X) would be an important step on the way to observing mixing with this 
technique. Recently, CLEO has observed a signal for — > K^tt~ (see Figure 2), and 
found R = B{D" K+TT-)/ B{D° K-tt+) ~ 0.8% [||. Unfortunately, without a preci- 
sion vertex detector, CLEO is unable to distinguish a potential mixing signal from DCSD. If 
the number of reconstructed charm decays can reach 10* around the year 2000, that would 
allow one to reach a new threshold of sensitivity to D^D^ mixing, and perhaps actually 
observe it. Therefore, it is time to take a detail look of all possible techniques for D^D'^ 
mixing search. 
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This paper |^ is organized as follows: in Section 2 there is a review of the experimental 
techniques which can be used to search for mixing, together with some thoughts on possible 
new techniques. In each case, the relevant phenomenology will be briefly presented. Section 
3 discusses the history, present status and future prospects of searching for mixing at different 
experiments. In Section 4, a comparison of the future prospects of the different experiments 
with different techniques, in the light of the CLEO II signal for JsT+tt", will be given. 

A brief summary is given in Section 5. Some detailed formulae and discussions (including 
possible CP violation effect) are provided in the appendices. 

2 The Techniques 

The techniques which can be used to search for mixing can be roughly divided into two 
classes: hadronic and semi-leptonic. Each method has advantages and limitations, which 
are described below. 

2.1 Hadronic method 

The hadronic method is to search for the decays 1)° K^tt^ {X). These decays 
can occur either through D^D^ mixing followed by Cabibbo favored decay _D" 
K^TT^ {X), or through DCSD D'^ K^tt^ (X). This means that the major complication 
for this method is the need to distinguish between DCSD and mixing |Q. The hadronic 
method can therefore be classified according to how DCSD and mixing are distinguished. In 
principle, there are at least three different ways to distinguish between DCSD and mixing 
candidates experimentally: (A) use the difference in the decay time-dependence; (B) use 
the possible difference in the resonant substructure between DCSD and mixing events in 
_D" — > iir+7r^7r",_ft'+7r^7r+7r^, etc. modes; (C) use the quantum statistics of the production 
and the decay processes. 

Method (A) requires that the be highly boosted and so that the decay time infor- 
mation can be measured. Method (B) requires knowledge of the resonant substructure of 
the DCSD decays, which is unfortunately something about which we have no idea at this 
time. Finally, method (C) requires that one use e^e~ annihilation in the charm threshold 
region. In the following, we will discuss these three methods in some detail. 

2.1.1 Method A -use the difFerence in the time-dependence of the decay 

This method is to measure the decay time of the _D" K^tt^ decay. Here the Z)" 
tagging is usually done by using the decay chain _D*+ -D"7r+ followed by — > K^tt^. 
The 7r+ from £)*+ has a soft momentum spectrum and is referred to as "the slow pion" . 
The idea is to search for the wrong sign Z)*+ decays, where the slow pion has the same 
charge as the kaon arising from the decay. This technique utilizes the following facts: 

(1) DCSD and mixing have different decay time-dependence, which will be described below. 

(2) The charge of the slow pion is correlated with the charm quantum number of the 
meson and thus can be used to tag whether a D'^ or l)" meson was produced in the decay 
D*+ D°Tr+ or D*' -D^tt". (3) The small Q value of the D*+ decay resuhs in a very 
good mass resolution in the mass difference AM = M{D*^) — AI{D^) — M(7r+) and allows 
a signal to obtained with very low background. (4) The right sign signal D^T^t 

^ This paper is essentially a revised version of Chapter 6 in and p^ . 
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Figure 1: The history of the quest for D^D'^ mixing. Note that the range in E691 result 
reflects the possible effects of interference between DCSD and mixing, and the CLEO II 
signal could be due to either mixing or DCSD, or a combination of the two. 
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Figure 2: The CLEO II signal for Z)° K'^n^. The D° mass for wrong sign events, (a) 
for events in the AM peak; (b) for events in the AM sidebands. The solid lines are the fits 
using the corresponding right sign mean and a in data. 
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followed by Z?" K 7r+ can be used to provide a model-independent normalization for the 
mixing measurement. 

A pure state generated at t = decays to the ii'+7r~ state either by D^D^ mixing 
or by DCSD, and the two amplitudes may interfere. The amplitude for a decays to 
if+7r~ relative to the amplitude for a decays to if ^7r+ is given by (see appendix A) 



A = ^Rmixing/2 t + VRdCSD 6^"^ (1) 

where 4) is an unknown phase, t is measured in units of average D° lifetime. Detailed 
discusion on the interference phase (f) can be found in Appendix A. Here Rdcsd — |pp 
where p is defined as: 

^~ Amp{D° ^ K+TT') ^ ^ 

denoting the relative strength of DCSD. We have also assumed a small mixing; namely, 
Sm,"f- <C 7+ or X, ?/ <C 1, and CP conservation. Detailed formulae and discussions (in- 
cluding possible CP violation effect) can be found in Appendix A. In the following, we will 
simply discuss the basic idea of how to distinguish DCSD and mixing with this technique. 

The first term, which is proportional to t, is due to mixing and the second term is due 
to DCSD. It is this unique attribute of the decay time-dependence of mixing which can be 
used to distinguish between DCSD and mixing. Now we have: 

1(1?" ^ K+TT-){t) OC (RdCSD + v/2RmixingRDCSD t coscj) + ^ Rmixi„gi^)e~* (3) 

Note that this form is different from what people usually use (but equivalent), see Appendix 
A. I prefer this form since it is not only more convenient for discussion here, but also much 
easier to be used to fit data. Define a = Rmixing/RDCSD, which describes the strength of 
mixing relative to DCSD. Equation ^ can then be rewritten as: 

1{D" K+'K-){t) cx Rdcsd(1 + %/2a tcoscf) + ^af)e-* (4) 

From this equation, one may read off the following properties p6[ |: (I) The mixing term 
peaks &i t = 2. (2) The interference term peaks at t = 1. (3) A small mixing signature can 
be enhanced by DCSD through interference (with coscj) ^ 0) at lower decay times, compared 
to the case without interference (with coscj) = 0). The ratio between the interference term 

and the mixing term, denoted £,{t), is given by ^(t) = \/ ^ cos(j>/t cx J i. So when a ^ 0, 



^ ^ oo. (4) Only for t > ^ ^\coscj>\ does the interference term become smaller than the 
mixing term. (5) I(to) ~ happens and only happens when coscj> = —1, and only at location 
— —. (6) One can obtain a very pure DCSD sample by cutting at low decay time. 



While Property (1) tells us that the mixing term does live at longer decay time. Prop- 
erty (3) tells us clearly that we should not ignore the interference term. In fact, that's the 
last thing one wants to ignore! (unless we know for sure cos cj> — 0). The commonly believed 
"annoying background", namely DCSD, could actually enhance the chance of seeing a very 
small mixing signal through the interference, compared to the case without the interference. 
In other words, the "annoying DCSD background" does not necessary inherently limit the 
hadronic method since the potentially small mixing signature could show up in the inter- 
ference term. For a very small mixing rate, almost all the mixing signature could show up 
in the interference term, not in the mixing term, as long as coscj) ^ 0. Property (2) tells 
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us at which location one expect to find the richest signature of a potential small mixing, 
which is where the interference term peaks: t ^ 1 (why should one keep worrying about 
long lived DCSD tails? let's hope for cos cj) 7^ first). Property (5) shows that destructive 
interference is not necessarily a bad thing. In fact, it could provide extra information. For 

example, if cos0 = — 1, then one should find I(<o) = at to = Figure |[ Note this 

unique attribute will become more interesting in method B, see Appendix B. This tells us 
that the destructive interference does not necessarily wipe out the signature of mixing. For 
the general case, interference will lead to very characteristic time distribution, as can be 
clearly seen in Figure Property (6) shows that we can study DCSD well without being 
confused by the possible mixing component. This will also become more important when 
we discuss method B. 

Therefore the signature of mixing is a deviation from a perfect exponential time distri- 
bution with the slope of 7+ |^ Our ability to observe this signature depends on the number 
of ^ K^TT~ events we will have. Right now this is limited by the rather poor statistics. 
Figures || or Figure ^ shows each term with a = 10% and cos0 = ±1 (with Rdcsd = !)■ 

It is worth to point out that the interference between mixing and DCSD also occurs 
in B^B'^ system. In this case, mixing is quite large and can be well measured while DCSD 
is small and unknown. The signature of the small DCSD would mostly show up in the 
interference term. But here we are not interested in measuring mixing nor measuring DCSD, 
what is interesting here is to measure CP violation. In Appendix C, We will use B^ — > D^tt~ 
as an example to show the basic idea. 

It is interesting to point out here that there is also a possibility, previously unrecog- 
nized, of using the Singly Cabibbo Suppressed Decays (SCSD), such as D° K^K^ ,tt^tt^ 
to study mixing ||2^ . This is because (assuming CP conservation) those decays occur only 
through the CP even eigenstate, which means the decay time distribution is a perfect ex- 
ponential with the slope of 71. Therefore, one can use those modes to measure 71. The 
mixing signature is not a deviation from a perfect exponential (again assuming CP con- 
servation), but rather a deviation of the slope from (71 -I- 72)/2. Since 7+ = (71 -I- 72)/2 
can be measured by using the D'^ —>■ K^tt^ decay time distribution, one can then derive 
y = 7_/7_|_ = (72 — 71)/ (71 +72)- Observation of a non-zero y would demonstrate mixing 
caused by the decay rate diflference (Rmixing — {x^ + y^)/2)- It is worth pointing out that in 
this case other CP even (odd) final states such as KsP^ can be also used to measure 

71(72)- In addition, there is no need to tag the , since we only need to determine the 
slope. Note that this method is only sensitive to mixing caused by the decay rate difference 
between the two eigen states, not to mixing caused by the mass difference x — Sm/"f-^- 
(Sm = m2 — mi). Right after this technique was introduced [ p6[ last summer, Fermilab 
fixed target experiments E791 started to apply this idea to their data j4^. The sensitivity 
of this method is discussed in Section 4.1. 

2.1.2 Method B - use difference in resonance substructure 

The idea of this new method |^ is to use the wrong sign decay D*^ D'^'^t followed 
by D'^ K^tt~tt'^, if+7r"7r+7r^, etc., and use the possible differences of the resonant 
substructure between mixing and DCSD to study mixing. There are good reasons to believe 
that the resonant substructure of DCSD decay is different from that of mixing (Cabibbo 
favored decay, CFD). We can use the D'~' K^tt^tt^ decay as an example. Detail discussion 

^One can use K^tt^ to study the acceptance function versus decay time. 
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Figure 3: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%. 
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Figure 4: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%, 
in log scale. 
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Figure 5: The decay time dependence of DCSD and mixing with maximal destructive 

interference cosip = —1.0. For different a = Rmixing/RDCSD vahics: from left to right, 
a = 100%, 50%, 10%, 5%,1%,0.5% (with Rdcsd = lO'^, this corresponds to Rmixing = 
10-2, 5 X 10-3, 10-3, 5 X 10-^, 10-^, 5 X 10-^). 
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Figure 6: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%. 
For different cos(/> values: from bottom to top, cos(f) = —1.0, —0.99, —0.96, —0.94, —0.80, 
—0.6,0.0,0.5,0.7, 1.0. The soUd line is the DCSD term, as a reference line. 
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about this method can be found in appendix B (including possible CP violation effect) , here 
we will just outline the basic idea. 

For CFD and DCSD, the true yield density n{p) at a point p in the Dalitz plot can be 
written as: 

n{p) cx l/i e'^^sfc + /2 e'^^-EWp+ip) + h e'^'BWK'-{p) + /4 e'^^ BWk.o{p)\^ (5) 

where fi are the relative amplitudes for each component and <f>i are the interference phases 
between each submode. is the S-wave three-body decay amplitude, which is assumed to 
be flat across the Dalitz plot. The various terms BW are Brcit-Wigner amplitudes for the 
_D" — > K*n and Kp sub-reactions, which describe the strong resonances and decay 

angular momentum conservation: BWr oc j^^. -_^MH-iTR/2 '^here Mr and are the mass 
and width of the My resonance {K* or p), and On is the helicity angle of the resonance. For 
CFD, fi and (pi have been measured by MARKIII |Q, E691 and are being measured 
by CLEO II. For DCSD, fi and 0i have not been measured. Note that in general 

f DCSD I r CFD I f DCSD , r CFD / ■ i ■\ fa\ 

h /fi r Jj ffj i^rj) (6) 

J DCSD I 1 CFD /-X 

This means that the resonant substructure (the true yield density n{p)) for DCSD is different 
from that of mixing. As both DCSD and mixing contribute to the wrong sign decay, the 
yield density for the wrong sign events riwip) will have a complicated form. Just like in 
method A, for very small mixing, the interference term between DCSD and mixing could 
be the most important one. 

Mathematically, the time-dependence of D° — > K'^tt~tt^ is the same as that of D° 
if+TT", the only difference is that now both Rdcsd and the interference phase ip (between 
DCSD and mixing) strongly depends on the location p on the Dalitz plot. As discussed in 
Appendix B, the time-dependence can be written in the form pq|: 



1 



«d(p) + V 2R.mixing nnip) nc{p) COS(/)(p) t + - nc{p) Rmixing t" 



e-^ (8) 



where riD {p) and nc [p) are the true yield density for DCSD and CFD respectively. Detailed 
discussion on the interference phase </> can be found in Appendix B. 

In principle, one can use the difference between the resonant substructure for DCSD and 
mixing events to distinguish mixing from DCSD. For instance, combined with method A, one 
can perform a multi-dimensional fit to the data by using the information on AM, M{D°), 
proper decay time t and the yield density on Dalitz plot riu, (p, t) . The extra information on 
the resonant substructure will, in principle, put a much better constraint on the amount of 
mixing. Of course, precise knowledge of the resonant substructure for DCSD is needed here 
and so far we do not know anything about it. Because of this, for current experiments this 
method is more likely to be a complication rather than a better method when one tries to 
apply method A to D° ^ K+tt-tt° (see ^ and ^) or D° K+tt-tt+tt-. In principle, 
however, one can use wrong sign samples at low decay time (which is almost pure DCSD) 
to study the resonant substructure of the DCSD decays. 

It is interesting to point out here (as discussed in detail in Appendix B) that the Dalitz 
plot changes its shape as the decay time "goes by" due to the interference effect. Note that 
the interference phase, unlike in the case of — > K^Tr~ , strongly depends on the location 
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on the Dalitz plot since there are contributions from the various Breit-Wigner amplitudes, 
which changes wildly across each resonance. One would expect that cos could have any 
value between [-1,1], depending on the location p. It is interesting to look at the locations on 
the Dalitz plot where cos (j>{p) = — 1 (maximal destructive mixing-DCSD interference). As 

pointed out in method A, Property (5) tells us that one should find I(io) = at = 
This means that the mixing-DCSD interference would dig a "hole" on the Dalitz plot at 
time to at that location. Since to oc \/Rdcsd(p) = \/ "-d (p) / nc (p) , the "holes" would show 
up earlier (in decay time) at locations where Rdcsd(p) is smaller. Imagine that someone 
watches the Dalitz plot as the decay time "goes by" , this person would expect to see "holes" 
moving from locations with cos (j){p) = — 1 and smaller Rdcsd(p) toward locations with 
cos (j){p) = — 1 and larger Rdcsd(p)- The existence of the "moving holes" on the Dalitz plot 
would be clear evidence for mixing. Once again we see the importance of the mixing-DCSD 
interference effect. 

In the near future, we should have a good understanding of DCSD decays and this 
method could become a feasible way to search for mixing (and CP violation) . 



2.1.3 Method C — use quantum statistics of the production and decay processes 

This method is to search for dual identical two-body hadronic decays in e^e" — s- 
^f" D^D^ , such as {K~ n^){K^ -k^) , as was first suggested by Yamamoto in his Ph.D 
thesis |5^. The idea is that when D^D^ pairs are generated in a state of odd orbital 
angular momentum (such as ^")^ the DCSD contribution to identical two-body pseudo- 
scalar- vector {D — > PV) and pseudo-scalar-pseudo-scalar [D PP) hadronic decays (such 
as (i4r~7r+)(i4r"7r+)) cancels out, leaving only the contribution of mixing The 
essence of Yamamoto's original calculation for the {K~ ■k'^){K~ tt^) case is given below. 

Let's define e^{t) = e^™'*-'^'*/^ ^ 2, 2) and e±{t) = {ei{t) ± e2{t))/2. A state 
that is purely \D'^) or |iD°) at time t = will evolve to \D{t)) or \D{t)) at time t, with 
\D{t)) = e+it)\D°) + e-(t)\D°) and \D{t)) = e-{t)\D°) + e+{t)\D°). In e+e^ ^ ^ 
D'-'D'^, the D^D^ pair is generated in the state D^D^ — D^D^ as the relative orbital angular 
momentum of the pair C — 1. Therefore, the time evolution of this state is given by 
\D{t)D{t')) - \D{t)D{t')), where t (f) is the time of decay of the D (D). Now the double- 
time amplitude Aw{t,t') that the left side decays to K~'k^ at t and the right side decays 
to K~^t:^ at t', giving a wrong sign event (K^'K^){K^Tr^), is given by: 

A^^{t,t') = (e+(t)e_(t') - e_(t)e+(t'))(a' - &') (9) 

where a = {K^Tr^\D'^) is the amplitude of the Cabibbo favored decay K^tt^, while 

b = {K'TT+\D°) is the ampHtude of DCSD D° K-tt+. Similarly, the double-time 
amplitude ^^(^,^0 for the right sign event [K^ 'k^){K^ n^) is given by: 

Ar{t,t') = ie+{t)e+{t') - e_(t)e_(i'))(a' - b') (10) 

One measures the wrong sign versus right sign ratio R, which is: 

N{K^TT+,K^TT+)+N{K+Tr~,K+TT~) _ JJ\An,{t,t')\^dtdt' 

~ N{K-7r+,K+n-)+N{K+n-,K-n+) ~ JJ\Arit,f)\^dtdt' ^ ' 

Note in taking the ratio, the amplitude term (a^ — b^) in Equations ^ and |l^ drops out. 
Thus, clearly R does not depend on whether b is zero (no DCSD) or finite (with DCSD). 
Integrating over all times, one then obtains R — (a;^ -f- y^)/2 = Rmixing, where x and y are 
defined as before. 
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This is probably the best way to separate DCSD and mixing. The exclusive nature of 
the production guarantees both low conibinatoric backgrounds and production kinematics 
essential for background rejection. This method requires one use e^e~ annihilation in the 
charm threshold region. Here the best final state is (if ~7r+)(i<r^7r+). In principle, one can 
also use final states like {K~ p^){K~ p'^) or (if *~7r+)(if *~7r+), etc., although again there 
are complications. For example, it is hard to differentiate experimentally {K^ p^){K~ p^) 
from (if ^(0+)(A'^7r+7r°), where DCSD can contribute. With high statistics, in principle, 
this method could be combined with method B. 

It has been pointed out that quantum statistics yield different correlations for the D^D'^ 
decays from e+e" D° D°, D°D°-f, D°D° Tr° Q. The well-defined coherent quantum 
states of the D^D^ can be, in principle, used to provide valuable cross checks on systematic 
uncertainties, and to extract x = Sm/^-^- and y = 7-/7+ (which requires running at different 
energies) if mixing is observed |^3| . 

2.2 Semi-leptonic method 

The semi-leptonic method is to search for Xl~v decays, where there is 

no DCSD involved. However, it usually (not always!) suffers from a large background due 
to the missing neutrino. In addition, the need to understand the large background often 
introduces model dependence. In the early days, the small size of fully reconstructed samples 
of exclusive D'^ hadronic decays and the lack of the decay time information made it difficult 
to constrain the D°D° mixing rate using the hadronic method, many experiments used semi- 
leptonic decays. The techniques that were used were similar — searching for like-sign fJ-'^JJ^ 
or /Lt-/i- pairs inp+N -> p+{p+p+)X |||,|35| and7r-Fe p+ n+ H, tt-W p+p+ 
These techniques rely on the assumptions on production mechanisms, and the accuracy of 
Monte Carlo simulations to determine the large conventional sources of background. 

There are other ways of using the semi-leptonic method. The best place to use the 
semi-leptonic method is probably in e+e^ annihilation near the charm threshold region. 
The idea is to search for e+e" ^ D°D° [K' 1+ v){K- 1+ v) or e+e" D- D*+ 

{K^'K~'K~){K^l~v)iT'^ The latter is probably the only place where the semi- 

leptonic method does not suffer from a large background. It should have a low background, 
as there is only one neutrino missing in the entire event, threshold kinematics constraints 
should provide clean signal. 

It has been pointed out that one can not claim a D^D^ mixing signal based on the 
semi-leptonic method alone (unless with the information on decay time of D^). Bigi [^Sj 
has pointed out that an observation of a signal on ~> l^X establishes only that a certain 
selection rule is violated in processes where the charm quantum number is changed, namely 
the rule ACharm = — AQi where Qi denotes leptonic charge. This violation can occur either 
through mixing (with the unique attribute of the decay time-dependence of mixing), 

or through new physics beyond the Standard Model (which could be independent of time). 
Nevertheless, one can always use this method to set upper limit for mixing. 

3 Mixing Searches at Different Experiments 

3.1 e+e running on *"(3770) MARK III, BES, Tau-charm factory 

The MARK HI collaboration was the first (though hopefully not the last) to use the 
e+e" ^" D^D^ technique. They reported three events consistent with |A5 = 2| 
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transitions |5^. One event is observed in the final state K~^Tr~ vs i^T+vr^Tr". The other two 
are reconstructed in the final states K^n^n'^ vs _fi'+7r~7r°, and a Dalitz plot analysis finds 
one to be consistent with K~ versus p'^ and the other consistent with K*'^t:^ versus 
K*°TT^ (note the I?" K-tt+tt^ decays are dominated by D° R- p+ and D° K*^tt^ 
channels). Using a maximum likelihood analysis, they interpreted the results for two limiting 
cases: a), if there is no DCSD in A'+tt^ and i4r+7r~7r°; then the events imply 

Rmixing = (1-2 ± 0.6)% or Rmixing > 0.4% at 90% C.L.; b). if there is no DCSD for 
D'^ K^Tr~ and no D^D^ mixing, and also at least one of the if +7r~7r°'s in each of those 
two events are non-resonant, then the results imply Rdcsd — r(Z?° — > K'^tt~tt'^)/T{D^ 
K'TT+TT^) = (7±4) tan'^6'c or Rdcsd > 1-9 tan'* 6*0 at 90% C.L.. This was a interesting 
result at that time, and had a strong influence on the subject. However, one cannot draw 
a firm conclusion about the existence of D^D^ mixing based on these events. There are 
at least two reasons: (1) The background study has to rely on Monte Carlo simulation of 
the PID (particle identification - Time-of-Flight) 0. As Gladding has pointed out: "These 
results must be considered preliminary because the calculation of the confidence level is 
sensitive to the tails of PID distribution for the background" [Q; (2) Assuming that the 
Monte Carlo background study is correct, and that the events are real, one still cannot claim 
the two events are due to mixing, for example, the non-resonant decays —> Kmr^ may 
contribute to one side of the pair in each of the events, in which DCSD can contribute. 

The MARK III puzzle can be solved at a r-charm factory, which is a high luminosity 
(10'^^CTO~^s~'^) e+e" storage ring operating at center-of-mass energies in the range 3-5 GeV. 
The perspectives for a D'^D'^ mixing search at a r-charm factory have been studied in some 
detail |5|]. I will outline here the most important parts. The best way to search for 
mixing at r-charm factory is probably to use e+e" — * ^" D^D^ [K^ tt^){K^ tt^). 
The sensitivity is not hard to estimate. Assuming a one year run with a luminosity 
of 10^'^cTO^^s^^, 5.8 X 10^ D^s would be produced from 5'". Therefore about 9 x lO'' 
(_ftr~7r+)(ii'+7r^) events would be produced. About 40% of them (3.6 x 10**) could be fully 
reconstructed. A study [ p7[ has shown that the potential dominant background comes from 
doubly misidentified {K^ tt^){K'^-k~), and if TOF resolution is 120 ps, this background 
could be kept to the level of one event or less. This means one could, in principle, set an 
upper limit at the 10"'' level. 

As mentioned Section 2.2, the best place to use the semi-leptonic method is probably at 
a r-charm factory. One good example is to search for e+e~ — > D^D*^ (_ft'+7r"7r+)(if+Z"i^)7r+. 
It is expected that this method can also have a sensitivity at the lO"** level. There are many 
other independent techniques that one can use for a mixing search at a r-charm factory. By 
combining several independent techniques (which require running at different energies), it 
was claimed that D'^D'^ mixing at the 10~^ level could be observable fs^ . 

There have been several schemes around the world for building a r-charm factory. If 
such a machine is built, it could be a good place to study mixing. The history of the r-charm 
factory can be fomid in reference |59|: one was proposed at SLAC in 1989 and one at Spain 
in 1993; one was discussed at Dubna in 1991, at IHEP (China), and at Argonne in 1994 
and at this workshop. It will be discussed again at IHEP (China) soon. Let us hope that 
we will have one in the not-too-distant future. 



^In principle, one can use kinematics to check whether the events are due to doubly 
misidentified D° K'TT^n^: if one inverts the ^ tt~ assignments and recalculates the 
_D" mass (let's call this Mflip), the MfHp will not be within mass peak if it is real, unless 
the K and tt momentum is close. Unfortunately, Mflip for all the three events are within 
mass peak. That's why one needs to totally rely on PID. 
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3.2 e+e running near T(4S^) ARGUS, CLEOII, CLEO III, B factory 



Without a precision vertex detector, CLEO II can only in effect measure the rate 
B{D^ Ktt) integrated over all times of a pure D° decaying to a final state Ktt. The ratio 
K—B{D'^ K^TT^) /B{D'^ K^TT^) is given by integrating equation || over all times (see 
appendix A) 

R. = R-mixing + R-DCSD + V 2RmixingRDCSD COScf). (12) 

CLEO finds || R = (0.77 ± 0.25 (stat.) ± 0.25 (sys.))%. This signal could mean one of 
two things: (1) mixing could be quite large, which would imply that mixing can be observed 
in the near future; (2) the signal is dominated by DCSD. The theoretical prediction for 
Rdcsd is about (2 — ?>)tan'^Oc ~ (0.6 — 0.9)% |6|, which is quite consistent with the 
measured value. It is, therefore, believed by many that the signal is due to DCSD, although 
it remains consistent with the current best experimental upper limits on mixing, which are 
(0.37-0.7)% @ and 0.56% [|8). 

CLEO has also tried to use hadronic method B, by searching for K'^tt~tt'^. The 

excellent photon detection at CLEO II allows one to study this mode with a sensitivity close 
to K~^TT~ mode. The main complication faced here is that (as discussed in method B) 

the resonant substructure is not necessarily the same for wrong sign and right sign decays. 
Because of this, as discussed in appendix B, the interpretation of R as Rmixing or Rdcsd will 
be complicated by the lack of knowledge of the details of the interference between submodes 
(and also the decay time information). Moreover, one has to worry about the detection 
efficiency across the Dalitz plot. Setting an upper limit for each submode is clearly very 
difficult. CLEO has set an upper limit ^ on the inclusive rate for K^it~it^ as 

R = B{D° K+TT-TT^) /B{D^ K-TT+TT^) < 0.68%. Note this upper limit includes the 
possible effects of the interference between the DCSD and mixing for each submode as well 
as the interference between submodes. 

This summer, CLEO will install a silicon vertex detector (SVX) with a longitudinal 
resolution on vertex separation around 75 /im. This will enable CLEO to measure the 
decay time of the , and reduce the random slow pion background (the resolution of the 
]^*+ _ £)0 mass difference is dominated by the angular resolution on the slow pion, which 
should be greatly improved by the use of the SVX). By the year 2000, with CLEO III (a 
symmetric B factory) and asymmetric B factories at SLAC and KEK, each should have 
thousands of — > K'^ K~ {X),'k'^'k~ {X) and a few hundred K'^Tr~ (and perhaps 

1)° i^+TT^TT*^, X+7r~7r+7r~ too) signal events with decay time information for one year of 
running. The typical decay length of 13° (£) is about a few hundred /im, and the resolution 
of the decay length (ac) is about 80 /im {C/ac ~ 3). The sensitivity to mixing at CLEO 
III and asymmetric B factories has not been carefully studied yet. A reasonable guess is 
that it could be as low as 10^*. If mixing rate is indeed as large as DCSD rate, it should 
be observed by then. 

3.3 Fixed target experiments 

A significant amount of our knowledge has been gained from Fermilab fixed target 
experiments, and in fact the current best upper limits on mixing have emerged from these 
experiments (E615, E691), and will come from their successors E687, E791 and E831 soon. 

The best upper limit using the semi-lcptonic method comes from the Fermilab exper- 
iment E615, which used a 255 GeV pion beam on a tungsten target. The technique is to 
search for the reaction TriV D'^D^ (K~ ^~^v)D^ {K~ ii^v){K~ fX^v), where only 
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the final state muons are detected (i.e. the signature is hke-sign ^jl^ i-i^ or /i pairs). 
Assuming ct(cc) ~ nuclear dependence, they obtained Rmixing < 0.56% Q. 

The best upper limit using the hadronic method by measuring the decay time infor- 
mation comes from E691, which is the first high statistics fixed target (photoproduction) 
experiment. In fact, E691 was the first experiment which used the decay time information 
(obtained from the excellent decay time resolution of their silicon detectors) to distinguish 
DCSD and mixing. The decay chains D*+ D°Tr+ followed by D° K+Tr', K+tt-tt+tt- 
were used. Their upper limits from the K^tt^ mode are Rmixing < (0.5 — 0.9)% 

and Rdcsd < (1-5 — 4.9)% , while the upper hmits from D° — > if+7r~7r+7r~ are Rmixing < 
(0.4 — 0.7)% and Rdcsd < (1-8 — 3.3)% . The ranges above reflect the possible effects of in- 
terference between DCSD and mixing with an unknown phase {(/>). Although the combined 
result gives Rmixing < (0.37 — 0.7)%, in principle, one cannot combine the results from the 
two modes since the interference phases are totally different, as discussed in appendix B. 
Note that in their analysis for A'+tt^tt+tt^, the resonant substructure in the Cabibbo 

favored and DCSD decays was ignored. As discussed in appendix B, in general, one can- 
not treat D'^ tt~ n'^ tt~ exactly the same way as — > A'+tt^ when one ignores the 
resonant substructure. Moreover, in principle, the detection efficiency vs decay time may 
not be studied reliably by using K^ir^ir^Tr^ as the resonant substructure could be 
different for CFD and DCSD. 

At the Charm 2000 workshop |63j, both E687 and E791 reported their preliminary 
result from part of their data. The best upper limits on mixing should come from these two 
experiments soon. Some most recent preliminary results can be found in p6t , progress has 
been made ^ on the measurement of the lifetime difference between _D" K~ii^ and 
-> K+K-{tt+tt-). 

4 Comparison of Different Experiments 

4.1 Hadronic method A 

This measurement requires: (1) excellent vertexing capabilities, at least good enough 
to see the interference structure; (2) low background around the primary vertex. The 
background level around the primary vertex could be an important issue as the interference 
term in Equation 3 does peak &t t = 1. In addition, low background around primary 
vertex means that one does not suffer much from random slow pion background and also 
one can measure the DCSD component at short decay times well. This also means that 
good acceptance at short decay times are very important. These are also important for 
understanding DCSDs at large decay times. The vertexing capabilities at e+e" experiments 
{C/a ~ 3) for CLEO III and asymmetric B factories at SLAC and KEK may be sufficient 
for a mixing search. The extra path length due to the Lorentz boost, together with the 
use of silicon detectors for high resolution position measurements, have given the fixed 
target experiments an advantage in vertex resolution (typically C/ a 8 — 10) over e+e~ 
experiments. One major disadvantage at fixed target experiments is the poor acceptance at 
short decay times. The low background around the primary vertex at e^e~ experiments is 
a certain advantage. It is worth pointing out here that at the e^e~ experiments (especially 
at an asymmetric B factory or Z factory) it maybe possible to use D*'^l~v, where 

the primary (Z?*+ decay) vertex can be determined by the /~ together with the slow pion 
coming from the Z?*+. In this case, the background level around the primary vertex is 
intrinsically very low ||2^ . 
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However, in the case of Z?" K^K^ , tt+tt^, etc., the requirement on the background 
level around the primary vertex is not so important. In this case, the mixing signature is 
not a deviation from a perfect exponential (again assuming CP conservation), but rather 
a deviation of the slope from (71 + 72)/2. It is worth pointing out that there are many 
advantages with this method. For example, one can use Cabibbo favored decay modes, 
such as K~Tr~^, to measure the average decay rate (71 + 72)/2 (which is almost a 

pure exponential, mixing followed by DCSD effect should be tiny, see Appendix A). This, 
along with other SCSD CP even (or odd) final states, would allow for valuable cross checks 
on systematics uncertainties. In addition, since we only need to determine the slope here, 
we do not need to tag the and do not have to use the events close to the primary 
vertex. The sensitivity of this method depends on how well we can determine the slope 
difference, which needs to be carefully studied. This is currently under study j4|, 
Roughly speaking, in the ideal case, the sensitivity to y would be ~ I/a/ZV, where N is the 
number of —> K~ ,tt^t:~ , etc. events, which means that the sensitivity to mixing 
caused by the decay rate difference (^ 2/^/2) would be close to 1/N. For example, a 
fixed-target experiment capable of producing ^ 10* reconstructed charm events could, in 
principle, lower the sensitivity to ~ 10~^ — 10~^ level for the term in Rmixing — (2;^+y^)/2. 
In reality, the sensitivity depends on many things and should be carefully studied. 

It is worth to point out that the current PDG experimental upper limit ( 90% CL) on 
the life time difference is only 

WrjO — TrjO I 

' ' = 2 y < 17% (13) 



This is based on the upper limit Rmixing = {x'^ + J/^)/2 < 0.37%, which is the E691 [|o) 
combined results on the K^n^ and D'^ A'+tt^tt+tt^ modes by assuming no 

interference between DCSD and mixing for both modes at the same time. 



4.2 Hadronic method B 

In the near future, we should be able to have a good understanding of DCSD Q in 
_D" — > Ar+7r~7r*', K~^tt~tt^tt~ , etc. modes, then method B will become a feasible way 

to study mixing and the sensitivity should be improved. Just like method A, this method 
requires very good vertexing capabilities and very low background around the primary 
vertex (this is even more important than in method A, since precise knowledge of DCSD 
is very important here). In addition, this method requires that the detection efficiency (for 
the mode being searched) across Dalitz plot be quite uniform (at least the detector should 
have good acceptance on the Dalitz plot at locations where DCSD and mixing resonant 
substructure are different). This is necessary so that detailed information on the resonant 
substructure can be obtained in every corner on the Dalitz plot. 

The excellent photon detection capabilities will allow e^e~ experiments to study the 
D° K+TT-n^ mode with very low background. From the CLEO II — > AT+tt tt" 



analysis [g7| 25 , the detection efficiency across the Dalitz plot will have some variations 
due to cuts needed to reduce background, however, it is still good enough to obtain detailed 
information on the resonant substructure. Future fixed target experiments may have a good 



^It may be possible that good understanding of DCSD can be reached by measuring the 
pattern of £)+ DCSD decays where the signature is not confused by a mixing component. 
It is worth pointing out that the £)+ DCSD decays can be studied very well at future fixed 
target experiments and B factories. 
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chance to study _D" — > K^tt^tt^tt^ mode, since the detection efficiency across Dahtz plot 
should be quite flat. The sensitivity that each experiment can reach by using this method 
depends on many things and need to be carefully studied in the future. 

4.3 Hadronic method C 

The sensitivity of this method depends crucially on the particle identification capa- 
bilities. Since the is at rest, the K and tt mesons will have the same momentum, so 
a doubly misidentified 13° — > A'^7r+ {K^ 7r^,7r+ K^) mimics a K^tt^ with 

almost the same D*^ mass. It is worth pointing out here that particle identification is not 
as crucial to method A as it is to this method (C), as far as this particular background is 
concerned. This is because in method A, the is highly boosted, and doubly misidentified 
K~^Tr~ decays will have a broad distribution in the D'^ mass spectrum around the 
D'^ mass peak; this background can be kinematically rejected with only a small reduction 
of the efficiency for the signal events. 

Once the sensitivity reaches 0(10^^), one may have to worry about other contributions, 
such as contributions from continuum background, contributions from e^e^ ^ 27 — )■ D^D^ 
which can produce C-even states where DCSD can contribute [ p2[ . 

4.4 Semi-leptonic method 

The semi-leptonic method usually suffers from large background (except at a t charm 
factory), the traditional method of looking for like sign /i^/i^ or /i^/i^ pairs is an example. 
New ideas are needed in order to improve the sensitivity significantly. Some promising 
techniques have been suggested by Morrison and others at the Charm 2000 workshop 
and have been discussed in the working group psf . 

The technique suggested by Morrison is very similar to that of the hadronic method: 
one uses the decay chain Z3°7r+, instead of looking for — > X+tt", one can search 

for — + K^l~v where there is no DCSD involved. Of course, due to the missing neutrino, 
this mode usually suffers from large background. However, for events in which the neutrino 
is very soft in rest frame, K^l~v is quite similar to K^tt~ kinematically. 

In this case, one has the same advantages as D*^ D^tt^ followed by K^tt~ has. In 

addition, as the neutrino is soft, the proper decay time of the can be reasonably estimated 
from K^l~ . The potential mixing signal therefore should show up as a term in the proper 
decay time distribution. To select the events with soft neutrino, one can require the K^l~ 
mass above 1.4 GeV. This requirement will keep about 50% of the total signal. One major 
background here is the random slow pion background, as the effective mass difference width 
is still much larger (a factor of 10) than Z3*+ D^tt^ followed by D° K^tt~ . In order to 
reduce this background, Morrison has suggested to look for a lepton with the correct charge 
sign in the other side of the charm decay. Another background is DCSD decay — ^ K^tt~ 
when the 7r~ fakes a l~ , however, this background will only populates at the higher end of 
the K^l~ mass spectrum where the neutrino energy is almost zero. This can be eliminated 
by cutting off that high end of the K^l~ mass. In principle, this idea can be used in a fixed 
target experiment as well as in a e^e~ experiment. The sensitivity of this method depends 
on the lepton fake rate (meson fakes as a lepton). One can find some detail discussions in 
Morrison's Charm2000 workshop summary paper ]6^ . 

Another technique, suggested by Freyberger at CLEO, is based on the technique which 
has been used by ALEPH, HRS and CLEO to extract the number of I3*+ D^T^t events. 
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The technique utiUzes the following facts: (1) Continuum production of cc events are jet like. 
(2) The jet axis, calculated by maximizing the observed momentum projected onto an axis, 
approximates the D*^ direction. (3) The D*^ -D"7r+ decay is a two-body process, and 
the small amount of energy available means that the 7r+ is very soft, having a transverse 
momentum p± relative to the D*~^ direction which cannot exceed 40 MeV/c. This low 
transverse momentum provides the _D*+ — > D^t:^ signature. 

The facts are used in the following way. The maximum momentum in the lab that 
the 7r+ can have perpendicular to the line of flight of the Z)*+ is 40MeV. One can define 
this quantity as pj_ = |p^|sm07r, where sm0^ is the angle between the D*+ and the tt^ in 
the lab frame, and is the magnitude of the 7r+ momentum. Hence, the 7r+ from 
will populate the low p±^ (or sinO.,^) region. The signal is enhanced if one plots p\ (or 
sin^OTr) instead of p±. One then looks for an lepton in the jet with the correct sign, namely, 
7r+Z+ right sign combination and tt+Z" wrong sign combination. The signal D*^ —> -D^tt^ 
followed by -D° K^l^v will peak in the low p^ (or sin^OT^) region for the wrong sign 
events. It is worth pointing out that one can look for a lepton in the other side of the event 
to reduce background. 

There are many kinds of background to this method one has to worry about. One of 
the major backgrounds is fake lepton background. For example, the decay chain _D*+ — > 
D^Ti^ {K~X)tt^ will also peak at the low p\^ (or sm^6'^) if the K~ is misidentified as a 
l~ . Another major background is probably the tt*^ dalitz and 7 conversions in — > Xn^ 
followed by tt" — > 7e+e~ or Xtt'^ followed by tt" 77 and then 7 e+e~. These 

two major backgrounds are at about 0.3% level in the current CLEO II data. Understand- 
ing these backgrounds is the major difficulty faced by this method. Although for CLEO 
III, things should improve, it is not clear what kind of sensitivity one can expect from this 
method for future experiments. Nevertheless, it is an interesting idea and worth investigat- 
ing. In fact, this technique is currently under study at CLEO [ |67| . 

5 Summary 

The search for D°D° mixing carries a large discovery potential for new physics since 
the D^D^ mixing rate is expected to be very small in the Standard Model. The past decade 
has seen significant experimental progress in sensitivity (from 20% down to 0.4%). Despite 
these 18 years of effort there is still much left to be done. 

As was discussed in the introduction, the observation of K^Tr~ is an important 

step on the way to observing a potentially small mixing signal by using this technique. With 
the observation of D° K+tt' signal at CLEO at the level R = 0.0077 ± 0.0025 (stat.) ± 
0.0025 (sys.), any experimenter can estimate the number of reconstructed Z?*^ — > K^n^ 
events their data sample will have in the future and what kind of sensitivity their experiment 
could have. Eighteen years after the search for D° K^n^ started, we have finally arrived 
at the point where we could take advantage of possible DCSD-mixing interference to make 
the mixing search easier. 

In light of the CLEO's _D" K^n~ signal, if the mixing rate is close to that of DCSD 
(above 10~^) , then it might be observed by the year 2000 with either the hadronic or the 
semi-lcptonic method, either at fixed target experiments, CLEO III, asymmetric B factories 
(at SLAC and KEK), or at a r-charm factory. If the mixing rate is indeed much smaller 
than DCSD, then the hadronic method may have a better chance over the semi-leptonic 
method. This is because the semi-leptonic method usually suffers from a large background 
due to the missing neutrino, while the hadronic method docs not. Moreover, the commonly 
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believed "annoying DCSD background" does not necessary inherently limit the hadronic 

method as the potentially small mixing signature could show up in the interference term. 
The design of future experiments should focus on improving the vertexing capabilities and 
reducing the background level around the primary vertex, in order to fully take advantage 
of having the possible DCSD and mixing interference. In addition, we have learned that 
the very complication due to the possible differences between the resonant substructure in 
many DCSD and mixing decay modes — > K^-k~{X) could, in principle, be turned to 
advantage by providing additional information once the substructure in DCSD is understood 
(the method B) and the sensitivity could be improved significantly this way. This means 
that understanding DCSD in D decays could be a very important step on the way to observe 
mixing. Experimenters and theorists should work hard on this. 

In the case of K'^tt~{X) and X^l~, we are only measuring Rmixing = 

{x^+y^)/2. Since many extensions of the Standard Model predict large x = 6m/^+, it is very 
important to measure x and y separately. Fortunately, SCSD can provide us information 
on y. This is due to the fact that decays such as ,ti^tt~ , occur only through 

definite CP eigenstate, and this fact can be used to measure the decay rate difference 
J/ = 7_/7+ = (72 — 7i)/(7i + 72) alone. Observation of a non-zero y would demonstrate 
mixing caused by the decay rate difference. This, together with the information on Rmixing 
obtained from other methods, we can in effect measure x. I should point out here that x and 
y are expected to be at the same level within the Standard Model, however we do not know 
for sure exactly at what level since theoretical calculations for the long distance contribution 
are still plagued by large uncertainties. Therefore, it is very important to measure y in order 
to understand the size of x within the Standard Model, so that when D^D^ mixing is finally 
observed experimentally, we will know whether we are seeing the Standard Model physics 
or new physics beyond the Standard Model. 

In this sense, it is best to think of the quest to observe mixing (new physics) as a 
program rather than a single effort. 
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A Appendix A — The Time Dependence of D'^ K^tt 

Note that Appendices A, B and C are from reference p5| . 

A.l The time-dependent effect 

A pure -D° state generated at i = could decay to K~^7t~ state either by — 
mixing or by DCSD, and the two amplitudes may interfere. Following the notation in psf , 
the time evolutions of \Di > and \D2 > are given by 



\D,{t) >= e^lA >, = e 



(* = 1,2) 



(14) 



with mi,^i {i — 1,2) being the masses and decay rates of the two CP eigenstates, The mass 
eigenstates \Di > and I-D2 > are given by 



where 



\Di >=p\D" > +q l^o > 
\D2 >=p\D° > -q > 

p _ 1 + e 
a " 1 -e 



(15) 
(16) 

(17) 



Under the phase convention \D^ >= CP\D'^ >, a state that is purely |D" > {\D" >) 
prepared by the strong interaction at t = will evolve to |Dpjjyg(t) > {\D°phys{t) >) 



phys(0 >- 2 



(ei+ 62)1-0° >+-(ei-e2)|i^o> 
P 



(18) 



l^°phys(t) > = 



(ei- 62)1-0" > +(ei+ 62)17^" > 



(19) 



Let us define a(/) = Amp{p" /), a{f) = Amp{D" /) with p{f)_^ a{f)/a{f); and 
a{f) = Amp{D^ ^ /), a(/) = Amp{D_^ /) with p(/)_= a{!)/a{f). Now the decay 
amplitude for states initially pure (D^) to decay to / (/) is given by (with / = K^tt~ 
and define |/>= CP|/ >) 



Ampd-OO (i) >^ /) 



«(/) 



^mp(|-DOphys(i) >^ /) - 2 °(/) 



+ 62) + - (ei - 62) 
P 



- (ei -62) + p{f) (ei +62) 



Therefore, we have 



/( \dXM >-> / 



— — |-| I (61-62) + -p(/) (61+62) I 



/( |i50phys(t) >->/) = (61 - 62) + ^Pif) (61 + 62) P 

A q p 



(20) 
(21) 

(22) 
(23) 



Note that neither | nor p{f) is an observable by itself. The flavor eigenstates and D'^ 
are defined by the strong interaction only, thus the relative phase between -O" and are 
undetermined: jZ)" >= e^*"CP|-D° > (we have set a — for convenience). The phase 
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transformation {D'^ >^ e''°'\D^ > leads to ^ ^ ^^'^ Pif) ~^ ^ *"p(/)- However, under 
the phase transformation, ^p{f) — > ^Pif) remains unchanged which means that ^pif) is 
an observable and we define r/ = ^pif) [rj = ^p{f) )■ 



We also have 



ei + 62 = 2 e" 



cosh(-(5m t H — t) 
2 4 



and 



ei 



62 = 2 e-™+*-^* sinh(-5mt + ^^7^) 

2 4 



(24) 



(25) 



where m+ = (toi + m2)/2, (5m = m2 ~ mi, (57 = 72 — 7i and 7_|_ = (71 + 72)/2 which is the 
average D° decay rate. With S, = ^{ iSm + ^S'f )t, we get 



I{ iDl^ysit) >- / ) = I^P I sinh(e) +7^ cosh(e) P e^-^^* 

/( 7^""phys(0 / ) = |a(/)P I^P I sinh(0 +^ cosh(e) \' e~^+'. 



(26) 
(27) 



Note for 5mt,S^t << 1 and for small |?7| ( in the case of / = JC+tt ), equation |26 
becomes 



/( \D"it) / ) = |a(/)|2 I i(*(5m + ^(57) t + P e-^+* 



1 



(28) 



^(l^physW /) = 

47]^ + ( {Sm f + 
and equation |2^ becomes 

/(|D%hy.(t) >^ /) = 



mr ,9, 



((57)^ 



)t^ + 2( SjRej] + 25m Imr;) i 



-7+t 



477^ + ( {5mf + )t2 + 2( (57Re77 + 2(5m Im^) t 



-7+t 



(29) 



(30) 



Note the difference between equation ^ and |30| is the indication of CP violation. These 
equations are in the form which people usually use. Next I will write the equations in a 
different form, which I believe is more convenient for discussion here. 

Now let us look at the same equation in a different way. For convenience, let us assume 
CP conservation so that we have |^| = 1 and |a(/)| = |a(/)| = \a\. Equation |2^ can be 
written in the form 



I{\D%M >-/) = |an ^^ + r;pe-* 



(31) 



where now the time t is measured in unit of average lifetime (1/7-1-), and 



(5m 



(5m — m2 — mi. 



V = 



7+ 



7± 



72 ± 71 
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Note again that we have assumed a small mixing; namely, Sm, 7_ ^ 7+ or a;, y ^ 1, which 

+ ?/^)/2. In addition, let's define Rdcsd = Note when 



means we have Rmixing 
If I ~ 1, Rdcsd = \v? 



which is the natural definition of Rdcsd (see equation 



Now we have 



/(|Z?;jhys(i) /) - lap I VRmixing/2 t + ^RdCSD ^ e" 



where 



(j) — Arg{ix + y) — Arg{r]) 
which is an unknown phase, and this gives: 

/(li^OhysW >^ /) - l«l' X 



Rdcsd + y/ 2RmixingRDCSD t COS(f>+ -Rmixing t 



(32) 
(33) 



(34) 



Note the unknown phase (f) ~ Argiix + y) — Argirj) depends not only on the mixing 
parameters, which are the sign and relative size of x {5m) and y ((^7), but also on the relative 
phase between DCSD D'^ K^tt^ and Cabibbo favored decay K^tt^ which is due 

to final state interaction. We can define ^mixing — Argiix + y) with (?!)dcsd = ^'"5(^) = 
Arg{p{f)). Note x = ^/2R~~ sin^mixing and y ^/2R~~ cos 



Arg{^^) 



'-'mixing ■ 



By 



using equation g4|, one can only measure Rdcsd, Rmixing and cos0 (up to a sign for (j)), but 
not (fiBCSD which depends on final state interaction and is unknown, neither (/>mixing which 
depends on the sign and relative size of x (Sm) and y ((^7). Note that since i] ^ rj under 
the phase transformation \D^ >^ e'"|D'^ >, 0dcsd = ^fgiv) is a physical parameter. In 
order to measure (/)dcsd and (/fmixing, one can use D° K'^K'^ etc. to measure y (^7) as 
discussed in method A. Together with Rmixing and cos 4> (up to a sign for ^) measured using 



equation |3J, one can thus in effect measure x or 5m (up to a sign) 
to a sign) and (/)dcsd(= ilfZ-mixingl ± \(l>\)- 



therefore 



-'mixing 



(up 



It has recently been argued by Browder and Pakvasa jT^I that it is possible that 0dcsd 
can be calculated. Although at present the calculations are purely phenomenological and 
are plagued by large uncertainties, it is worth investigating since it would be very helpful if 
0DCSD can be calculate theoretically in a reliable way. 

As discussed above, the difference between equation ^ and |3^ is the indication of CP 
violation. Equation 29 and ^ can be written in the form: 



I{\D%M >- /) = I«(/)P l-P X 



P 



Rdcsd + \/ 2Rmixing Rdcsd t cos(j)+ -Rmixing 



(35) 



and 



J(|i50p,,,(t) >^/) = |a(/)p l^p X 



Rdcsd + y 2RmlxlngRDCSD t COS0+ -Rmixing i 



(36) 



where Rdcsd — |?7p and = Arg{ix + y) — Arg{fj). Note 4> is different from (p because 
Arg{fj) = Arg{^) + Arg{p{f)) where Arg{^) — ~Arg{^). For convenience, let us ignore the 
CP violation in the decay amplitude. Thus the difference between the interference terms 
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(that is, the interference phase (j) and 0) in equation ^ and |3^ would be the indication of 
CP violation. 

It has been recently argued [^8[ ^ that CP violation effect should not be ignored 
here, since new physics which can introduce large mixing rate may often involve significant 
CP violation. In the SU(3) limit, the quark model and factorization gives Arg{p{f)) — 0. 
In this special case, and also assuming that x » y (large mixing caused by new physics) 
which means (^mixing = 7r/2, the interference term now simply becomes ± xsin(2ipM) t 
where ^ = e~^"^*'. Under those specific assumptions, the interferece term is odd with 
respect to CP, thus CP violating. This is the case discussed recently in |69]. In general, 



due to final state interaction, we have |50| Arg{p{f)) ^ 0. In addition, in the Standard 
Model, X and y are expected to be at the same level. But we do not know exactly at what 
level since theoretical predictions are still plagued by large uncertainties. For instance, if 
Rmixing ~ 10"'' in the Standard Model Q, then y should be somewhere between 10""^ to 
10~^. Not to mention the current upper limit on y is at 10~^ level (this is why measurement 
on y by using — > K^K^ etc is so important, as discussed before). Therefore, unless new 



physics introduces large x far above 10""^, one should use equation ^5| and 36 or equation |29| 
and ^ to fit data, instead of assuming x » y (or (/)mixing = ^/2)- It is worth to point out 
here that we are looking for a tiny effect, the possible contribution due to the decay rate 
difference y should not be ignored unless x » y is experimentally confirmed. 

As discussed in method A, a small mixing signature could be enhanced by DCSD 
through interference at lower decay times, compared to the case without interference. This 
really depends on the actual value of cos 0. If nature is unkind, the sign and relative size 
of X (dm) and y (^7), and 0dcsd caused by final state interaction are such that cos0 = 
COs((^jnixing (/•dcsd) — 0, then we will not have this advantage by using this technique. In 
this case, we may have to use K^tt'tt^ (method B) where coscj) strongly depends 

on the location on the Dalitz plot, and cannot be zero everywhere on the Dalitz plot. In 
general, K^ir^ir'^ may provide much more information on mixing and CP violation 

than — !■ K^TT^ does. We will discuss this point in more detail in appendix B. 

It is worth to take a look at the time dependence of the right sign decay K~^Tr~ 
here. The decay amplitude for states initially pure to decay to / is given by: 



Ampi\D%i,ys{t) >^f)^\ Hf) 



P 

- p(/)(ei - 62) + (ei + 62) 



(37) 



Therefore we have 

/( |i5%hys(i) >^ / ) - h (ei - 62) + (ei + 62) \'' (38) 



which is 

/( D%bys{t) >^ / ) = |a(/)P \v sinh(e) + cosh(e) P 6"^+*. (39) 
Similarly, we have 

/( D°hy.(i) >->/) = l«(/)P 1^ sinh(e) + cosh(e) \' e-^^K (40) 
Note for 5mt, 5^t << 1, equations ^ and ^ become 

/( \D%^ys{t) >^ / ) = |a(/)|2 I 1 rj {tSm + ^Sj) t + 1 \' 6"^+' (41) 

and 

/( |Z?^hysW / ) - l«(/)P I I V + U-f) t + 1 p e-^+' (42) 
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Assuming CP conservation, Equation ^ can also be written in the form: 



/(|C"ph.vs(t) >^ /) = \aV I VR-DCSDRmixing/2 te^* + IY e 



where 



(f) = Arg{ix + y) + Arg{ri) 
which is an unknown phase, and this gives: 

Ii\D%i.ys{t) >^ f) = lap X 

1 + v/2R mixing Rdcsd t cos(j) + -R 

mixing Rdcsd t 



(43) 
(44) 



(45) 



From equation 45 one can see that the mixing term is further suppressed by Rdcsd due 
to the fact that K^tt~ can occur through mixing followed by DCSD K^Tr~. 

The interference term is similar to that of equation |3^, but the interference phase (/) is 
different from the (p in equation 34, Since the first term is the dominated one, equation k3 



is essentially the same as /(|D''phys(i) >— /) 



A. 2 The time-integrated effect 

It is interesting to take a look at the time integrated effect. Let us define 

jrniDXsit) >- f)dt 



R = 



li\D%^ys{t) >^f )dt 



and we have 



R = Rdcsd + V 2RmixingRDCSD cos (j) + Rmixing- 

In the special case when | cos0| — 1, equation^ can be written in the form 
( Rdcsd — R )^ + ( Rmixing — R )^ = R^. 



(46) 
(47) 

(48) 



If one can measure R and Rdcsd precisely, then Rmixing can only have two possible 
values Rmixing = R± \/ 2Rdcsd (R — Rdcsd /2) where the sign ± corresponds to cos (j) — ^1. 
For example, if one finds R = Rdcsd = Ro, then Rmixing can be either or 2Ro. 



In general, we do not know cos 0, so we cannot determine R; 



mixing 



this way. How- 
ever, there are still useful information in the time integrated effect. For instance, one can 
compare R{t < 0.2), which is measured by using events that have decay times t < 0.2, 
to the R for all decay times. At short decay times where mixing has not yet fully devel- 
oped, the wrong sign events will be almost pure DCSD (see Section 6.2.1), which means 
we have Rdcsd — R(i < 0.2). Equation |4^ tells us that R will be different from Rdcsd 
by 2RmixingRDCSD COS (j) + Rmixing- Compared to Rdcsd, the fraction of the difference 

is A/2Rmixing/RDCSD COS -|- Rmixing/RDCSD- For Rmixing ~ 10"'', RdCSD 7.7 X 10"^ 

and I cos(/)| = 1, this corresponds to a ± ~ 16% change, which could be measurable with 
reasonable amount of data. Note here that one can only measure AR = R — Rdcsd- Ob- 
servation of a non-zero AR would demonstrate the existence of mixing. In addition, a 
precise measurement on Rdcsd using wrong sign events at short decay times would help the 
time-dependent analysis using equation For instance, it could be helpful to fixed target 
experiment if Rdcsd has been already measured very well at other experiments. This is 
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because that at fixed target experiment, the detection efficiency at short decay times is very 
low, leading to a poor constraint on the large DCSD term and making the time-dependent 
analysis very difficult. 

As the detection efficiency cancels out in the wrong sign vs. right sign ratio ^, the 
method discussed above would be useful when one does not know the detection efficiency 
vs. decay time therefore cannot perform a reliable time-dependent analysis. However, 
in the case of K'^tt~, one can always use the right sign events K~tt~^ to 

study the detection efficiency. Nevertheless, this should still be a useful cross-check for the 
time-dependent analysis. For example, one can use the fitted cos0 value from the time- 
dependent analysis to determine Rmixing in the time-integrated analysis. This would be a 
valuable cross-check for the measured value of Rmixing for the time-dependent analysis. 

If one wants to measure CP violation, one could then compare equation for 
K~^Tr~ and K'n'^ separately. Note that in this case one does not need to measure 

Rdcsd which is good for experiments without decay time information at all. However, since 
one has to measure K^tt^ and K^n^ separately the detection efficiency would 

not drop out completely. This means it could be quite difficult since we are looking for a 
tiny effect here. 



B Appendix B — The Time Dependence of I?° K^n 7r° 

B.l The time-dependent Dalitz analysis 

In the case of K^n^n'^, the situation is very similar to that of K^n^. 

For a given location p on the I?" K~^tt~tt'^ Dalitz plot, the decay amplitude for states 
initially pure D'^ to decay to / is given by (with / = K~^tt~tt^) 



Amp{\Dl^^,{t) >^ f){v) = \-a{f,p) 



p(/,p)(ei + 62) + ^ (ei - 62) 



(49) 



where a(/,p) = Amp{D" ^ f){p), a{tp) = Amp{D° f){p) and with p) = a{f,p)/a{f,p). 
Note that mathematically, equation ^ is exactly the same as equation |2^. Therefore, for a 
given location p on the Dalitz plot, we have 

/( |i^°hys(i) f)ip) = Hf,p)\' sinh(0 + ^{p) cosh(0 P e-^+* (50) 

where again C = ^( iSm + ) t, and 7y(p) = gP{f,p) which now depends on the location 
p on the Dalitz plot. Note in general, |77(p)| ^ 1 does not hold. For instance, |?7(p)| could 
be large at locations on the Dalitz plot where DCSD amplitude is large due to constructive 
interference among DCSD submodes, while CFD amplitude is small due to destructive in- 
terference among CFD submodes. Nevertheless, as long as \'q{p)\ <C min{{^)~^ , i^^)^^), 
for Smt, S-ft << 1 we still have 

just like in the case of D'^ K^tt^ , this can be written in the form 



/( \D%^M >- / )ip) = Hf.P)\' I 7;i^Sm+-Sj) t + r,{p) f e"^-* (51) 



( |Z?;^hys(0 >^ f)iP) = Hp)\^ I v/Rmi-ng/2 ^ + VRdCSD (p) C^'^^P^ ^ e"* (52) 



^This is because we will combine and D* decays, and we have assumed CP 

conservation. 
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where we have assumed CP conservation, with |a(/,p)p — |a(/,p)p = |a(p)p, andRDcsD(p) 
|77(p)p as the definition of Rdcsd(p) at location p on the Dalitz plot. Now we have 

H Hi^ysit) >-/) (p) = Kp)P X 



RdCSd(p) + y 2RmixingR.DCSD(p) t COS 0(p) + -Rmixing t 

Note again the unknown phase = Arg{ix + y) — Arg{ri{p)) depends not only on the 
mixing parameters, which are the sign and relative size of x (Sm) and y (Sj), but also on the 
relative phase between DCSD D'^ K^tt^tt^ and Cabibbo favored decay — > i4r~7r+7r° 
at location p on the Dalitz plot. As in the case of K~^Tr~ , we can define 

^mixing = Arg{ix + y) 
p 

0dcsd(p) = Arg{r]{p)) = Arg{-) + Arg{ p{f,p) ) 

q 

4>{P) = ^mixing - 0DCSD {p) 



Let = |a(/,p)p = nD(p) which is the true yield density on the Dalitz plot for 

direct DCSD decay D'^ K'^tt^tt'^, and |a(p)p = nc{p) which is the true yield density for 
the direct CFD decay K^tt^t:'^, equation p3 can be written in the form 



/( \Dl^^M >-> / ) (P) - I v/R.ni-ng/2 t + e^'^(P) p e-\ (54) 

which gives 

J( |I?°hys(i) >^ / ) (P) = 



1 



2Rmixmg nD{p) nc{p) COS(/)(p) t+ - nc{p) Rmixing t" 



(55) 



where nc{p) and nD{p) can be written in the form 

nc{p) « I /i e'^^Aaf, + ^ e'^-BWp+{p) + h e'^^BWK>-{p) + h e'^'BWuAp) ?, (56) 



and 



31 e'«iA36+52 e'^-EWp 



93 e 



'BW. 



K'- + 



94: e 



'BW, 



(57) 



where fi {gi) are the relative amplitudes for each component for CFD (DCSD) and 
are the relative interference phases between each submode for CFD (DCSD). A^b is the 
S-wave three-body decay amplitude, which is assumed to be flat across the Dalitz plot. 
The various BW{p) terms are Breit-Wigner amplitudes for the K*-k and Kp sub-reactions, 
which describe the strong resonances and decay angular momentum conservation: BWr oc 
M --MR~iTR/2 ■^here Mu and F/j are the mass and width of the Afy resonance {K* or p), 
and On is the helicity angle of the resonance. 

The idea of the time-dependent Dalitz analysis is that, in principle, one could perform 
a multi-dimensional fit to the data by using the information on the mass difference, 
mass, proper decay time t and the yield density on the wrong sign Dalitz plot. Of course, 
in order to do this, one needs a very large amount of clean data - a formidable experimental 
challenge! Lest any reader's despair, we would like to remind the reader that CP violation 
was first discovered 31 years ago in Ki^ — > tt+tt^ with only 45 events with S/N ~ 1, and 
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now we are talking about 10 x 10^ A'l tt+tt events at KTeV with a expected background 
level at 0.03% 0! 

As we have pointed out, the extra information on the resonant substructure will, in 
principle, put a much better constraint on the amount of mixing (compare to that of — > 
if+TT^). A Monte Carlo study is currently underway and will be described elsewhere. Here 
we will just make some qualitative remarks: 

(a) The amplitudes /,; and relative interference phases in equation thus nc{p), 
can be measured very well for the Cabibbo favored decay D° K~tt^tt°. Note that only 
relative phases are measurable, so out of the four phases there is one unknown (denoted i^o)- 

(b) As pointed out in method A, in principle, one can use wrong sign samples at short 
decay time (which is almost pure DCSD) to study the resonant substructure of the DCSD 
decays. This means the amplitudes gi and relative interference phases 9i in equation ^t\ , 
thus nuip), can also be measured for the DCSD K^tt^tt^ . Note again that only 
relative phases are measurable, so out of the four phases there is one unknown (denoted ^o)- 

(c) Now let us look at ri{p) — ^p{f,p) where p{f,p) — a{f,p)/a{f,p). Assuming CP 
conservation, we have r]{p) = -b{p)/a{p) where |-| = 1 and the phase of - is unphysical 
as it depends on phase convention (we can set it to zero for convenience). Once we know 
fi{gi) and 4>i{9i) as can be measured in (a) and (b), then (/)dcsd(p) = ^''ff('?(p)) will only 
depend on (j)Q — 9q for a given location p on Dalitz plot. The sign and size of (j)o — do is due 
to final state interaction and is unknown. At different locations, the phase 0dcsd(j') will be 
different as it strongly depends on the location due to the phases of various Breit-Wigner 
amplitudes (BWr oc m- --'\iH-iT r/2 )- Nevertheless, in principle, for given CFD and DCSD 
resonant substructure, there is only one unknown phase which is <j)Q — Oq, and (po — Oq does 
not depend on location p on the Dalitz plot. 

(d) The phase ^mixing = Arg{ix + y) depends on the sign and relative size of x and y. 
For given x and y, ^mixing is fixed. 

(e) Therefore, for given CFD and DCSD resonant substructure and mixing parameters 
{x and y), the mixing-DCSD interference phase 4>{p) — Arg{ix + y) — Arg{rj{p)) will only 
depend on (pa — 9o at location p. For a given 0o — ^o, one would expect that cos (f>{p) could 
have any value between [—1, 1], depending on the location p. Therefore, one can perform a 
Monte Carlo study by changing (pQ — Oq and look at the wrong sign Dalitz plot at certain 
decay time t. This Monte Carlo study is currently underway. 

(f) It is interesting to look at the locations on the Dalitz plot where cos (j){p) = — 1 
(maximal destructive mixing-DCSD interference). As we have discussed in method A, at 
decay time to = a/ 2Rdcsd (p)/Rmixing , we have | v^Rmixing/S t + VRdcsd(p) e"^(P) | = 0. 
This means the mixing-DCSD interference would dig a "hole" on the Dalitz plot at time 
to at that location. Since to oc \/Rdcsd(p) = \/ n-u (j>) / nc (p) , the "holes" would show 
up earlier (in decay time) at locations where Rdcsd(?') is smaller. Imagine that someone 
watches the Dalitz plot as the decay time "goes by" , this person would expect to see "holes" 
moving from locations with cos (j){p) = — 1 and smaller Rdcsd(p) toward locations with 
cos (j){p) — —1 and larger Rdcsd(p)- The existence of the "moving holes" on the Dalitz plot 
would be clear evidence for mixing. 

(g) In practice, one can only expect to find the "holes" at earlier decay times, since at 
longer decay times there are not many events left. As already pointed out in (f), the "holes" 
showing up at early decay times have small Rdcsd(p)- In general, small Rdcsd(p) occurs at 
locations where the DCSD amplitude is small due to destructive interference among DCSD 
submodes while the CFD amplitude is large due to constructive interference among CFD 
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submodes. This means that it could be difhcuh to see the "holes" since there is not many 
wrong sign events to begin with at these locations, unless one has a very large amount of 
clean data. 

(h) There are also extra information at locations where cos (p{p) ^ —1. Therefore, the 
best way would be to perform a multi-dimensional fit to data by using the information on 
the mass difference, mass, proper decay time t and the yield density on the wrong sign 
Dalitz plot. 

In the discussions above, we have assumed CP conservation. As we discussed in ap- 
pendix A, CP violation effect could be important here. It is interesting to point out here that 
the extra information on the resonant substructure could also provide additional informa- 
tion on CP violation when one compares the difference between the resonant substructure of 
D'^ A'+7r^7r° and D'^ AT^tt+tt". As discussed in (c), assuming CP conservation, once 
we know fi{gi) and 4>i{9i) as can be measured in (a) and (b), then 0dcsd(p) = Arg{ri{p)) 
will only depend on (f>Q — 9o for a given location p on Dalitz plot. With CP violation, then 
Arg{r]{p)) will also depends on Arg{^). This is for Z?" K+ir^ir^. For 1?° Ar^Tr+vr", 
Arg{fi{p)) will depends on Arg(^) = ~Arg{^). This means the interference phase of 
— ^ Ar^7r+7r° will be shifted by 2 Arg{^ ) everywhere on the Dalitz plot, compared 
to that of D'^ — > Ar+vr^vr*'. This would lead to a totally different resonant substructure. For 
example, similar to what has been discussed in (f), if the person sees a "hole" at location 
p and decay time t on the K~^tt~tt'^ Dalitz plot, but could not find a "hole" on the 

Z)0 K~TT^n^ Dalitz plot at the same location p and decay time t, that would be indica- 
tion of CP violation. This will be studied by Monte Carlo soon. In any case, with enough 
data in the future, one should fit — > K^tt~tt^ and D'^ AT^tt+tt" separately, or study 
the difference between the two. 

As in appendix A, it is also worth to take a look at the time dependence of the right 
sign decay AT+vr^Tr" here. Assuming CP conservation, just as Equation we have 



I{\D%^y,{t) >^/) = |a(p)|2 X 



1 + -\/2RinixingRDCSD(p) t COS [p) + -R 

mixing R.dcsd(p) t 



e-* (58) 



where 



4> (p) = Arg{ix -l- y) -|- Arg{'r]{p)) (59) 



Unlike in the case of AT+tt" where Rdcsd ^ 0.77% which is very small, in the 

case of A'+7r~7r° Rdcsd (p) depends on location p and, in principle, could be quite 

large at certain locations. If mixing is indeed large, then the interference term could be 
important at locations where Rdcsd (p) is also large. For instance, if Rmixing ~ 10"'^ and 
for locations where Rdcsd (p) ^ 10, then the interference term ^ 2RmixingR.DCSD (j>) would 
be at 10% level. If mixing is indeed small, then the second and third terms can be ignored, 
and equation m simply becomes /(|I?°phys(i) /) = |a(p)P 

B.2 The time-integrated Dalitz analysis 

One alternative would be to study the time integrated Dalitz plot. Now equation |5^ 
becomes: 

^(l^phy. /) (P) « Hp)\' X 
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Rdcsd(p) + V2R mixing R.dcsd(p) cos(/>(p) + R 

n 



(60) 



As we have already discussed, at short decay times, mixing has not yet fuUy developed, the 
resonant substructure will be almost due to DCSD alone. Therefore, in principle, one can 
measure Rdcsd(p) or nD{p) this way. Note again that (j){p) = Arg(ix + y) — Arg{r]{p)), 
which can be written in the form 

</)(p) = + (61) 

where ^{p) is the part which depends on the location p on the Dalitz plot, and </> is the part 
which does not. 

Once we know the resonant substructure for CFD and DCSD, we know Rdcsd(p) and 
$(p). There are only two parameters in equation |60| which are unknown: (1) Rmixing which 
is what we are trying to measure; (2) which depends on Arg{ix + y) and (j)Q — 9q which is 
the overall relative phase between DCSD and CFD. Therefore, in principle, one could also 
perform a multi-dimensional fit to data by using the information on the mass difference, D'^ 
mass, and the yield density on the wrong sign Dalitz plot to measure these two unknown 
parameters. However, in practice, since we expect that the interference term is much larger 
than the mixing term, it may not be easy to separate Rmixing and the phase (jj effects since 
they enter the interference term as ^/Rmixing cos(<l'(p) + if)). 

One could also compare the Dalitz plot for decay times t < 0.2 (which is almost due to 
pure DCSD) to the Dalitz plot for all decay times (or for t > 0.2). At longer decay times, 
for given Rmixing, the resonant substructure would be modified according to cos0(p). Since 
cos (t>(j)) only depends on the location p, but not on the decay time t, the interference effect 



will not be integrated out. Equation 60 tells us that at location p the wrong sign resonant 



substructure will be changed by approximately 2RmixingRDCSD (p) cos(/)(p), compared to 
Rdcsd(p) measured at short decay times. The fraction of change is 2Rmixing/RDCSD {p) cos 
For Rmixing ^ 10^**, Rdcsd (p) ^ 10^^ and |cos0(p)| = 1, this corresponds to ± ~ 10% 
change, which could be measurable with reasonable amount of data. Note that since 
Rdcsd(p) is different at different location on Dalitz plot, the fraction of change will also be 
different at different location for given mixing rate. A clear difference between the Dalitz 
plot for decay times t < 0.2 (which is almost due to pure DCSD) to the Dalitz plot for all 
decay times (or for t > 0.2) would be clear evidence for the signature of mixing, although the 
information on the size of Rmixing could be somewhat washed out by the unknown phases. 

Note that the time- integrated Dalitz plot for K^tt^tt^ and — > K^tt^tt^ can 

also be used to study CP violation as discussed in Section 7.1. This is simply because the 
interference effect on the Dalitz plot will not be integrated out. 



B.3 The time-dependent Non-Dalitz analysis 

Another possible alternative would be to study the time-dependent effect after inte- 
grating over the whole Dalitz plot, or put it another way, one could ignore the resonant 
substructure and treat it simply exactly as K~^tt~ time-dependent analysis. In fact, 

this is what E691 did to i^+7r~7r+7r~ (E791 and E687 are doing the same thing now). 

With limited data size, this is a natural thing to try. However one must be very careful 
here since one cannot treat multi-body decays exactly as 1?° K^n^ . We will discuss this 
point below. 

After integrating equation [55| over the whole Dalitz plot, the decay rate /( |£'phys(^) 
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/ ) is now proportional to 



nc{p) da X 



/ nY){p)da 
J nc{p)da 

- ^ J V n^jp) nc{p) cos0(p) da 
' / nc{p)da 



Now one can define 



7?.DCSD 



/ 7iD(p)rfa 



(62) 



(63) 



/ nc{p)da 

which is what one could measure if there is no mixing (or use events at short decay times). 
However, in general, the interference term is not proportional to y/TZucSD- Therefore, 
one cannot treat D° if+7r~7r° (thus also I?" — *■ K~^tt~tt^tt~) exactly the same way as 
Z)° K~^TT~. Note that since the interference term in K^tt" and K~^tt~tt^ 

are totally different, one cannot combine the two results to set upper limit, as E691 did 
to D*^ — > K^Tr~ and D'^ — > K^tt^tt^tt^ by assuming the interference phase to be zero in 
both cases. 



In principle, if the following inequality holds ^ 



nu{p)da / nc{p)da> i a/ nD(p) nc{p) |cos0(p)| da 



(64) 



one can force 



/ V "d(p) nc{p) cos (j){p) da 



J nc{p)da 

With this definition, equation pq can be written in the form 



V'^DCSD cos^' 



ncip) da 



T^DCSD + v^2Rniixing T^DCSD t COS ^ + 2 ^mixing 



(65) 



(66) 



the 



Mathematically, equation 66 is similar to that of D'^ K^tt , i.e. equation But 
"interference" phase 5* has a totally different meaning since cos 5* strongly depends on 



the actual value of / y/ nD{p) nc{p) cos (/)(p) da. For example, in the extreme case where 
for any given value of ^y nY){p) nc{p), if cos 0(p) has an equal chance of being negative and 
positive (with the same size), then the interference effect could be completely integrated 
out, leading to cos^ = 0. In general, since cos4>{p) can be positive and negative, the 
interference effect could be greatly reduced. Therefore, this technique does not necessary 
have the advantage as in the D'~' — > K~^tt~ time-dependent analysis. 



B.4 The time-integrated Non-Dalitz analysis 

Another possible alternative would be to study the effect after integrating over decay 
time and also integrating over the whole Dalitz plot, this is in fact what has been done at 



CLEO l27l |2 



^ This inequality should hold in general. Intuitively, inequality ^ is similar to inequality 
S"=i SLi yj ^ ( SiLi V^iUi cos(/)j )^, which is trivial to prove. 
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Integrating equation |6^ over all decay times, and normalizing to the corresponding 
right sign samples, we get 



_ JnB{p)da /— — / a/ nujp) nc(p) cos0(p) da 

j nc(p)da V B J nc[p)da 

Note the information on mixing is essentially lost due to the lack of precise knowledge on 
riD (p) and cos (f>{p) . 



Again if the inequality 64 holds, equation can be written in the form 



n = 7^DcsD + v/2R mixing T^DCSD COS + Rmixing ■ (68) 

B.5 Summary of Appendix B 

Although it is true that more-than-two body hadronic decays, such as Z?" 
K~^Tr~Tr'^, are very complicated due to the possible difference in the resonant substructure of 
the DCSD and CFD (mixing) decays, this unique attribute could, in principle, provide addi- 
tional information which could allow one to distinguish DCSD and mixing, and to study CP 
violation in the future. The best way would be to perform a time-dependent Dalitz analysis, 
that is, one can perform a multi-dimensional fit to the data by using the information on 
AM, M{D^), proper decay time t, and the yield density on Dalitz plot nw{p, t). The extra 
information on the resonant substructure will, in principle, put a much better constraint on 
the amount of mixing and CP violation. For experiments without decay time information, 
one could perform a time-integrated Dalitz analysis. This could be a feasible way to observe 
mixing, but may not be a good way to measure or constrain mixing rate. With limited data 
size, one could in principle, perform a time-dependent non-Dalitz analysis, as one would 
do to D'^ — > K^TT" . However, care must be taken here since one cannot treat multi-body 
decays the same as — > X+7r~. In addition, since the extra information on the resonant 
substructure has been integrated out, the time-dependent non-Dalitz analysis would not 
have much advantage. Without decay time information and with limited data size, one 
could only perform a time-integrated non-Dalitz analysis as has been done in . 

As one can see, precise knowledge of the resonant substructure for DCSD is important 
to this technique. This means that understanding DCSD in D decays could be a very 
important step on the way to observing mixing and CP violation using this technique. In 
principle, one can use the wrong sign sample at very low decay times (which is almost pure 
DCSD) to study the resonant substructure of the DCSD decays. It is also worth pointing out 
that it may be possible that a good understanding of DCSD could be reached by measuring 
the pattern of -D"*" DCSD decays where the signature is not confused by a potential mixing 
component. In the near future, we should have a good understanding of DCSD decays and 
this method could become a feasible way to search for D^D^ mixing and CP violation. 



C Appendix C The Interference Between Mixing and DCSD in i?° Case 

As we have discussed, the interference between D^D^ mixing and DCSD can be used 
to study D^D'^ mixing and CP violation. It is interesting to point out that the interference 
between mixing and DCSD also occurs in B'^B^ system. We will use 5° D^tt~ as an 
example (for convenience, let's also call it DCSD, as it is doubly Cabibbo suppressed relative 
to i?° D^TT^). In this case, mixing is quite large and can be well measured while DCSD is 
small and unknown. The small signature of DCSD would mostly show up in the interference 
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term. But here we are not interested in measuring mixing nor measuring DCSD, what is 
interesting here is to measure CP violation. 



Equation ^ and ^ also apply to B'^ D^ti . We have (with / = D+tt ) 



Hf)\ 



1 |9|2 



sinh(^) + r\ cosh(^) 



'7+t 



(69) 



where C = 5 (idm + i (57 ) t. We have defined 77 ~ ^p{f), with p{f) = a{f)/a{f) and 
a(/) = AmpiB'^ ^ /), d{f) - Amp{BO ^ /). 

We also have 

/( B%hys(0 >->/) = \a(M 1^1' I sinh(e) +fy cosh(0 f 6"^+* 



(70) 



where we have defined 77 = j;p{f), with p{f) — a{f)/a{f), a{f) — Amp{B° — > /) and 
d{f) - Amp{BO /). 

In the case of B^^B^d mixing, 5^ <C 5m, which means we have ^ — ^Sm. Note that 
sinh(ia;) = i sin[x) and cosh(ia::) = cos(x), so that equation ^ and ^ becomes 

2 |9|2 



and 



J( |S° (t) >^/) = |a(/)n^|2 X 



|ryp + (1 - 77^) sin^(-(5TO t) + /m(ry) sin((5m t) 



J(i3%hy.(i) >-/) = K/)P X 

\fj\^ + (1 - sin^(i(5TO t) + /m(77) sin((5TO t) 



-7+t 



(71) 



-7+t 



(72) 



For simplicity, let us first neglect the phase difference between B^ D^tt and B^ 
_D+7r^ caused by final state interaction. Define 77 = \ri\ e** with Rdcsd — and Rdcsd = 



|77p, assuming |2| ~ 1, we have rj = rj* — \r]\ e 



/( / ) = l«(/)P X 

Rdcsd + (1 - Rdcsd) sin^(ia; t) + V Rdcsd sin$ sin(a;t) e"*, (73) 



and 



/( B%^ys{t) >->/) = |a(/)|2 X 

Rdcsd + (1 - Rdcsd) sin2(ia; t) - ^R-dcsd sin* sin(a;t) e"*, (74) 

where now the time t is measured in unit of average B^ lifetime (1/74-), and x = 

We do not know Rdcsd for B^ — > D+tt", but one would expect very roughly Rdcsd 
tan^Oc ~ 0.3%. Note again we have neglected the phase difference between B^ D^tt^ 
and B^ — + D^tt~ caused by final state interaction. Assuming sin$ = 1.0 (corresponds to a 
maximal CP violation), with x = 0.71 and Rdcsd ~ 0.3%, and also set |a(/)p = |fl(/)P — Ij 
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Figure ^ shows the decay time dependence of equation and together with each term 
alone. One can see that there is a sizeable interference between B^B'^ mixing and small 
DCSD. More importantly, there is a clear difference between /( \B'^Yiys{^) D'^tt~ ) and 
J( B^p\^ys{t) >— ^ _D^7r+ ). Note that the mixing term peaks at i ~ 2, where mixing is fully 
developed and the maximal asymmetry occurs. Figure ^ shows the same thing, but with 
sin<I> = 0.5. 

It is interesting to take a look at the time-integrated decay rate: 



r( iSpVs >-/) = l«(/)px 



R-DCSD + (1 ^ R-dcsd) 



2(1 + x2) 1 + 



VRdcsd sin$ 



(75) 



and 



r( I^Vys >^/) = K/)l' X 



R-DCSD + (1 ^ Rdcsd): 



DCSD 



sin$ 



'2(l + a;2) 1. 

Therefore the time- integrated asymmetry is (with Rdcsd ^ 2(i+x^) ) 
< Asym >— 

r( \B%,ys >^ D+TT- ) - r( IBOphys >^ D-n+ ) 



(76) 



^ 2 

X 



phys 



= >- 



D+TT- 

VRdcsd sin$. 



r( |S%hys >^ D-t:^ 



(77) 



With X = 0.71 and Rdcsd ^ 0.3%, we have < Asym >= 15% x sin<i> which is quite large. 
Similar conclusion can be made on the lepton-tagged D^tt^ on T(4S'). 

Of course, we have neglected the phase difference between i?'^ — > D+tt" and J5° 
Z?+7r^ caused by final state interaction. In general, however, one cannot ignore the phase 
difference and things would be more complicated. 

After I thought about this, I found that this mode has been discussed in the past. For 
people who are interested in this mode, relevant references can be found in to 

It is also worth to point out that there could be a possibility of using multi-body decays 
such as -B" Z?+7r~7r'' to study CP violation, similar to the case of I?'^ K^'it~tt^ where 
one could use the extra information on the resonant substructure. 
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Figure 7: The decay time dependence of /( l-Bphys(^) >^ D~^n ) (sin$ = 1.0) and 
/( BOphys(i) >^ D-7r+ ) (sin$ = -1.0) 
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Figure 8: The decay time dependence of /( l-Bphys(^) >^ D~^n ) (sin$ = 0.5) and 
/( BOphys(i) >^ D-7r+ ) (sin$ = -0.5) 
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